Predicted field-induced hexatic structure in an ionomer membrane 
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Coarse-grained molecular-dynamics simulations were used to study the morphological changes 
induced in a Nafion®-like ionomer by the imposition of a strong electric field. We observe the 
formation of novel structures aligned along the direction of the applied field. The polar head 
groups of the ionomer side chains aggregate into clusters, which then form rod-like formations 
which assemble into a hexatic array aligned with the direction of the field. Occasionally these lines 
of sulfonates and protons form a helical structure. Upon removal of the electric field, the hexatic 
array of rod- like structures persists, and has a lower calculated free energy than the original isotropic 
morphology. 

PACS numbers: 



The transport of protons through an ionomer mem- 
brane is at the heart of the operation of a polymer- 
electrolyte fuel cell. Any procedure that facilitates the 
transport would thus be of significant value in improv- 
ing the usefulness of this technology. In this Letter we 
describe simulations that indicate that application of a 
strong electric field could change the morphology of an 
ionomer membrane in such a way as to enhance proton 
transport appreciably. 

The study of proton transport in aqueous ionomers has 
received considerable attention [H, 0, Q . For perfluoro- 
sulfonated ionomer membranes such as Nation® Q the 
structure is characterized by a long, hydrophobic, flu- 
orinated main chain, and short side chains terminating 
in hydrophilic sulfonate anion groups. The dissimilarity 
between the hydrophobic backbone and the hydrophilic 
side chain terminations creates a phase-segregated mor- 
phology with sharp interfaces between the domains. The 
sulfonic-acid functional groups aggregate to form a hy- 
drophilic domain that is hydrated upon absorption of 
water. It is within this continuous domain that ionic 
conductivity occurs when protons dissociate from their 
anionic counterions and combine with water to become 
mobile hydronium ions. 

Beginning with a seminal paper by Eisenberg [5[, a 
number of authors have suggested morphological ge- 
ometries for this microphase separation. These include 
inverted globular micelles interconnected by a chan- 
nel structure of cylindrical micelles [6], and lamellar, 
sandwich-like structures. Other recent approaches are 
based on micelle-channel @, H[ and inverted micelle mod- 
els @. There is, however, no complete consensus as to 
the structure of the percolating networks of sulfonate or 
water channels. 

We have performed simulations of almost-dry Nation® - 
like ionomers with the goal of determining the mor- 
phological changes that would be produced by strong 



externally- applied electric fields. We find that the side 
chains self-organize into cylindrical clusters having their 
axes parallel to the applied field, and that these clus- 
ters form a hexatic array in the plane perpendicular to 
the field direction. At very strong fields each of the rod- 
like clusters has an ordered inner core consisting of dis- 
tinct wires of sulfonate head groups and their attendant 
protons. The stable cylindrical cluster structures that 
emerged included an unusual spiral arrangement of three 
sulfonic wires with four protonic chains. The structural 
changes resulting from the poling process appear irre- 
versible. 




FIG. 1: (Color online) This snapshot shows the ionomer be- 
fore application of a field. Large/small spheres represent the 
end- group oxygen/sulfur atoms of the side chains. The colors 
of the spheres represent altitudes, with blue at the base and 
red at the top of the simulation box. The size of all structural 
elements is schematic rather than space-filling. 
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Model - We employ a united-atom representation for 
Nafion 0, El EI El, where the CF 2 and CF 3 groups of 
the backbone and side chains, and the sulfur atom S and 
oxygen group O3 of sulfonate head groups are modeled 
as Lennard- Jones (LJ) monomers with a universal diam- 
eter o~ = 0.35 nm. We assume the number A of water 
molecules per sulfonate is small (A < 5), and that a pro- 
ton detached from a sulfonate captures a water molecule 
and becomes a hydronium ion. 

The force-field parameters for the ionomer are chosen 
to agree in most instances with the Nafion model of Pad- 
dison [3] and are given in our recent paper [l5[. All the 
partial charges on the ether oxygen, carbon and fluorine 
atoms of the side chain, and on the fluorocarbon groups 
of the backbone skeleton are set to zero. Electrostatic 
charge is located on the sulfur atom (charge +e) and 
oxygen group (charge — 2e), such that the total charge 
of the sulfonate head group SO3 is — e. The hydrogen 
ion H + is considered as part of a hydronium complex of 
charge +e. For the low water contents (A < 5) considered 
in this work, the remaining water molecules are mostly in 
an immobile state near the pore walls. Therefore one can 
effectively integrate out their degrees of freedom by intro- 
ducing a distance-dependent dielectric permittivity [l5[ 
and an effective proton diameter a. Aside from these two 
approximations, our approach is equivalent to the consid- 
eration of a membrane with A = 1, when all its sulfonic 



groups are dissociated forming hydronium ions 




FIG. 2: (Color online) View from above of the simulation 
system for applied fields of (a) 1.2 x 10 6 kV/m and (b) E = 
8.2 x 10 6 kV/m. A gain, a blue color represents a low altitude 
and a red color a high altitude. Protons are shown as black 
dots in diagram (b). 

The total potential energy of the membrane is taken 
to be 

u(^=J2 u t+lL u e+Y, u ™+Y, u ^ k -n\) (1) 

i j m k,l 

where (ri, tn) are the three-dimensional position 

vectors of the N particles in the system. In Eq. (pQ), 
i runs over all bonds, j runs over all bond angles, m 
runs over all torsional angles, and fc, I run over all non- 
bonded (Lennard- Jones and Coulomb) force center pairs 
in the system. For lightly humidified membrane materi- 
als, the protons are mostly energetically bound to their 



host sulfonate groups [H,[lH. The resulting SO3H effec- 
tive dipoles have zero net charge, and thus experience no 
net force in a uniform applied external field, but are sub- 
ject to torques. Interactions between dipoles are screened 
with a dipolar screening length R s [201 ]. 

In direct-current experiments, the electrostatic poten- 
tial U GX t(r) inside the polymer is considered to be linear 
[211 ]. In non-equilibrium molecular dynamic computer 
simulations the current flow is usually induced by ap- 
plying an external electrostatic field E j22[. The scalar 
protonic conductivity x ls defined from the linear rela- 
tion j = xE between the applied field E and the induced 
current density j in the membrane. When the external 
field is applied along the £-axis, j = ^qiVi yZ /V where 
Vi iZ is the axial component of the velocity of the z-th pro- 
ton, i = l,..,A/g, with N s the total number of protons 
and V the sample volume. The conductivity can also be 
related to the diffusivity D through the Nernst-Einstein 
relation x=Nse 2 D /VksT in equilibrium molecular dy- 
namics simulations. However, this method is less accu- 
rate for phase separated systems, where confinement ef- 
fects, such as proton trapping in the water-sulfonate clus- 
ters, strongly affects the overall diffusion constant D. 




FIG. 3: (Color online) (a) Depiction of only the backbone 
polymer shows the hexatic array of proton channels formed 
by the sulfonates. Circles and solid lines are guides for the 
eye. (b) Remanent structure in a membrane after removal of 
the applied field. 

At low A we expect surface hopping to be the lead- 
ing contribution to the current through the membrane. 
In this case the protons move in the electrostatic energy 
landscape of the sulfonate ions. Proton transfer is ex- 
pected to occur at sulfonate-sulfonate separations 0.7- 
0.8 nm 12311 with an energy barrier less than 2.1 kcal/mol 
I23I I24 |25| . We neglect the contributions to the 
current arising from Grotthuss and en-masse proton dif- 
fusion, which are only significant in highly humidified 
membranes (A > 5), and whose evaluation requires full 
simulations with explicit water molecules. 

Simulation results - Molecular dynamic simulations 
were performed for coarse-grained membranes kept at 
constant volume V =(11 nm) 3 and constant temperature 
T=300 K. Each side chain contains 2 hydrophilic and 7 
hydrophobic monomers, and there were 14 hydrophobic 
backbone monomers between adjacent sidechains. The 
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FIG. 4: (Color online) The xy plane projection (a) and a 
perspective view (b) of a single rod-like cluster with a five- wire 
structure from Fig. [2jb) . Small spheres attached to oxygens 
represent sulfurs and netted small balls are protons in the 
right figure. 




FIG. 5: (Color online) The xy plane projection (a) and a 
perspective view (b) of a single rod-like cluster with a spiral 
three- wire structure from Fig. EJb). Note that there are four 
protonic wires for a triplet of sulfonic wires. 
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FIG. 6: (Color online) Protonic conductivity % of the mem- 
brane as a function of applied field E. Line without /with 
circles depicts increasing/decreasing field. 

system temperature was controlled by a Langevin ther- 
mostat with a friction coefficient 7 = 0.1 and Gaussian 
noise 6ksTj. The equations of motion were integrated 
using the velocity Verlet algorithm with a time step of 0.2 
fs. Periodic boundary conditions and Lekner summation 
[2(^ of long-range electrostatic interactions were used. 
We temporarily detached the side chains from the 




FIG. 7: (Color online) A short-time average of the 3D density 
distribution of protons for material in zero field (a) and with 
£ = 8.2 x 10 6 kV/m (b). 



backbone skeleton [23 and cut the backbone into 14- 
monomer segments [28| for initial equilibrating runs of 50 
ps, after which the polymer was reassembled and equi- 
librated with another 50 ps run. Statistically averaged 
quantities were gathered during the next 5-10 ns. A snap- 
shot of a simulated system in the absence of an applied 
field is given in Fig. [TJ The ionomer has undergone a 
partial phase separation on the nanometer scale into a 
polymer phase consisting of backbone with pendant side 
chains and a hydrophilic phase formed by loosely con- 
nected clusters of the sulfonate head groups, which are 
mostly composed of compact multiplets of about 12 sul- 
fonates. 

Application of a strong electric field of the order of 10 9 
V/m gives rise to the formation of the new morphology 
illustrated in Fig. O The isotropic system of clusters in 
Fig. [T] is now replaced by a system of cylinders aligned 
along the z-axis. The sulfonates form the walls of these 
cylinders, which have a diameter of ~5<r (« 1.8 nm) and 
a separation of roughly 8-10cr. As the field strength is 
increased, these cylinders form a hexatic array, as shown 
in Fig. [3](a). At the same time, the internal structure 
of the cylinders self-organizes into an array of distinc- 
tive 'wire-like' chains of sulfonates, or to be more pre- 
cise, chains of oxygens and sulfurs, and parallel chains of 
protons. Though the number of sulfonic wires per cylin- 
der varies between 1 and 5, nearly half of the structures 
have a three-wire structure. All of the other structures 
observed consist of an equal number of sulfonate and pro- 
ton wires, an example for a five- wire structure being given 
in Fig. [U The mean SO3 -SO3 group spacing along the 
chains is about 0.65-0.80 nm, which is close to the sepa- 
ration distance expected for Nafion-like membranes [231 ]. 
The cylinders composed of three sulfonic wires accom- 
modate four protonic wires, as shown in Fig. [5] for the 
special case of a spiral cylindrical cluster. The appar- 
ent stability of this structure may be an artifact of the 
periodic boundary conditions. 

A remarkable result of the simulations was that this 
field-induced morphology remains after the removal of 
the applied field. Fig. [3](b) shows the ordered structure 
that persists after the field has been removed from the 
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system shown in Fig. [2](b). In order to understand this 
phenomenon we calculated and compared the free ener- 
gies A of the original isotropic structure and the struc- 
ture with remanent order. The method of thermody- 
namic integration tells us that the difference AA be- 
tween a system with and without Coulomb interactions is 
equal to J^~^ (dH(£)/d£) when the coupling parame- 
ter £ characterises the strength of Coulomb interactions 
U l J between head groups i and j. The configurational 
part of a parametrized Hamiltonian was chosen as 

17(C, r) = K + U{f) + (£ - 1) Ei>i ^ Th e calculated 
values of AA/NksT — —9.3 for the untreated membrane 
and AA/NksT = —9.7 for the treated membrane indi- 
cate that the latter has a lower free energy. The morpho- 
logical changes induced by the strong external field thus 
appear to be irreversible. 

Proton conduction is strongly enhanced by the forma- 
tion of the sulfonate cylinders along the current direc- 
tion, as indicated in Fig.[6l which shows that the effective 
conductivity x(E) increases rapidly with field strength. 
Because the effect of high fields is irreversible, there is 
hysteresis in the current-field plot. Both the small- field 
conductivity x = 9.5xl0 -3 S/cm of the isotropic mem- 
brane and the value ofx = 2xl0 -2 S/cm for the poled 
membrane are much smaller than the reported experi- 
mental values x — 0-1 S/cm for humidified Nation, but 
much larger than the theoretically predicted surface hop- 
ping conductivity x ~ 10 -5 S/cm [29|. 

The 3D density distribution of protons, plotted in 
Fig. shows a disordered distribution of protonic clouds 
inside the isotropic material. In the presence of a strong 
applied field, however, the protonic clusters connect to 
each other and form large clusters. These clusters elon- 
gate along the field direction, eventually forming rod-like 
channels. 

We thus conclude from our simulations that applica- 
tion of a strong applied electric field to a Nafion®-like 
ionomer may induce the formation of rod-like structures 
of sulfonate groups, which assemble into a hexatic array 
aligned with the direction of the field. This novel mor- 
phology persists after the removal of the field as a conse- 
quence of the lowering of the free energy associated with 
this change in structure. The induced agglomeration of 
hydrophilic head groups into long rods appears to reduce 
the percolation threshold for ion conductance [Io[. Sim- 
ulations also indicate that strong correlations between 
the sulfonate head groups and protons inside the rods 
create wire-like structures of charged particles. The ge- 
ometry of these structures is very rich and includes stable 
spiral-like three-wire clusters. These effects, if confirmed 
experimentally, may have an impact on the industrial 
applications of Nafion-like membranes. 
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